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Abstract 

In this paper we study the spectral gap for a family of interacting particles systems 
t;;j- . on [—N, N], proving that it is of the order A^~^. The system arises as a natural model 

CN ' for current reservoirs and Fick's law. 
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1 Introduction 



c^ ■ In this paper we study an interacting particle system whose state space is {0,1}[~^'^1. 

The dynamics is a Markov process with generator L = Lq + L^, Lq the generator of the 
stirring process (see (12. ip below), Lf, the generator of a birth-death process whose events 
are localized in a neighborhood of the end-points, see (13.11) . 

In particular we consider the case when around A^ there are only births while around 
— A^ there are only deaths. The system is then "unbalanced" and in the stationary measure 
fi% there is a non-zero steady current of particles flowing from right to left. The system is 
designed to model the Fick's law which relates the current to the density gradient. 

In statistical mechanics non-equilibrium is not as well understood as equilibrium, hence 
the interest from a physical viewpoint to look at systems which are stationary yet in non- 
equilibrium: in our case the stationary process is in fact non-reversible and the stationary 
measure fi% not Gibbsian. 

There is a huge literature on stationary non-equilibrium measures, in particular on 
their large deviations as they are related to "out of equilibrium thermodynamics" (see for 
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instance [H 121 [3l [8] ) . Here we study the spectral gap in the stationary process and the way 
it depends on the system size A^. Spectral gaps have been well studied in the reversible or 
Gibbsian set-up, both for stirring and for more general interacting particle systems. (See 
for instance [in]-) The techniques used in those situations however do not seem to apply 
to our non-equilibrium model. We shall rather rely on stochastic inequalities and coupling 
methods, thus reducing the problem to that of bounding the extinction time of the set 
of discrepancies between two coupled evolutions. The case of a single discrepancy can be 
regarded as an environment dependent random walk with death rate which also depends 
on the environment. Its extinction time has been studied in [7] and as we shall see here it 
is closely related to the spectral gap in our model. 

The main part of this paper refers to the case of "current reservoirs" (where also L^ is 
multplied by a factor 1/A^). Much simpler is the case when Lf, fixes the different densities 
at the boundaries, whose analysis is carried out sketchily in the next section mainly as an 
introduction. 



2 Density reservoirs 

We consider in this section the Markov process on {0, 1}[~^'^1 with generator L = Lq +L', 
where, denoting by r] the elements of {0, 1}[~^'^1, 

7V-1 

W(r/):=2 E [/(^^^'^^'^) - /(^)] (2-1) 

x=-N 

with r]^^'^'^^\x) = r](x + 1), 7^(^'^+i)(x + 1) = r](x) and ry(^'^+i) = r] elsewhere; 
L'fiv) = P+[/(^^^'^^)-/(^)] + (l-P+)[/(^^-'^^)-/(^)] 

+ P-[/(^(+'"^^) - fiv)] + (1 - P-)[/(^(-'-^^) - fiv)] 

where 1 > /0+ > P- > and ri~^'^{x) = 1, rj^'^iii) = viv)^ V ¥" ^j analogously r]~''^{x) = 0, 
ry-.^(y) = r]{y), y ^ x. 

The process corresponding to L' alone leaves unchanged the occupations at |x| < N 
while the equilibrium probabilities of occupation at ±A^ are equal to p±. Since p+ > P-, 
this creates a density gradient and the full process with generator L = Lq + L' describes 
the particles flux determined by the density gradient. The process is uniformly Doblin, 
in particular there is a unique stationary measure fi% and convergence to equilibrium is 
exponentially fast. The averages yU^[?7(a:)] describe a linear density profile in agreement 
with Pick's law. Fluctuations in equilibrium are well characterized, [11], and the large 
deviations as well, [S]. 

Denote by /ijv the initial distribution and by fJ^NSt the distribution at time t (i.e. the 
law at time t of the process with generator L starting from fiN)- Then, since the process 



is uniformly Doblin, for any positive integer A^ there are strictly positive constants cn and 
h^ so that 

WhnSi - IJ^%\\ < catc"''^* for any ^^ and t > (2.2) 

where for any signed measure A on {0, 1}[~^'^1 

l|A||=$^|A(r^)|. (2.3) 

We shall prove that: 

Theorem 1. There are c and b> independent of N so that for any initial measure fi^- 

\\^^NSt-^i%\\<cNe-'''"' (2.4) 

Proof. Let 

A'^v = {^Z = (^^'U^'^) e ({0,l}x{0,l})[-^'^l:r/^(a;) := r]^^\x) - r]^^\x) > 0,Va;}, (2.5) 
and, for / : X^ — )■ M, 

N-l 

^ofiv) := 2 E [/(!/^^'^^'^) - /(!Z)] 

x=-N 

+ P-[/(!/^+'"^^) - fiv)] + (1 - P-)[/(!Z^~'"^^) - fiv)] 

where r/(+'^)(x) = (1, 1), r/(-'^)(x) = (0, 0), x = ±A^. 

It is easy to see that Cq and C define Markov generators on A'^r. Moreover, when 
acting on functions that depend on only one of the two entries, t]*-^-* or r]^'^\ of r/, we see 
that Cq + C coincide with L, and so it defines a coupling between the processes with 
generator L starting from two comparable configurations r/^^^ and t]*^^-* {r]^^\x) > rj^'^^x) for 
all x), showing that the L-evolution is attractive in the sense of [2] (i.e. preserves order). 
In particular we may take r]^^^ = 1 and rj^^^ = the configurations that are identically 
1 and, respectively, 0. Moreover Cq leaves unchanged the number of discrepancies which 
instead may decrease under the action of C. Write P for the law of the process starting 
from 77*^^^ = 1 and t]*^^-* = and call 7r(x,t) = P[r]^(x,t) = 1]. We then have, recalling that 
7r(x, 0) = 1 for all x, 

7r(x, t) = l- I (ps{x, N)7r{N, t) + p,(x, -A^)7r(-A^, t)) ds 



where Ps{x, y) is the probabihty under the stirring process (with only one particle) of going 
from X to y in a time s; this is the same as the probability of a simple random walk whose 
jumps outside [— A^, A^] are suppressed. Hence 

7r(a:,t) = E4e-^*W] 

where E^ is the expectation of the process with transition probabilities ps{x,y) and T*(t) 
is the time spent at {— A^, A^} during [0,t]. We shall prove that 

^4e-^*W] <ce-^^"'* (2.6) 

which will then imply 

TV 

Y, P[T/^(x,t) = l] < A^ce-'^"'* 

x=-N 

and so (12. 4p . because fiNSt and fi% are squeezed in between the laws of the marginal of 
the coupled process. 
Proof of fl2i|) . 

By an iterative argument it is enough to show that 

sup E,[e-^]<p<l, t:=T*{N^) 

xe[-N,N] 

But 

inf PJt> 1] > 5>0 (2.7) 

as the probability of reaching {— A^, A^} by time A^^ — 1 is bounded from below uniformly 
in the starting point and the probability of not moving for a unit time interval is also 
bounded away from 0. By (12. 7p 

E,[e-^] = E,[e-^;T<l]+E,[e-^;r>l] 

< 1 - P4t > 1] + P4r > l]e-^ < 1-6(1- e~') 

n 



3 Main result 

In this paper we study the process with generator L = Lq + Lb, Lq as in (12. ip . Li, = 
Lf)^^ + L(,^„ describes births and deaths near the boundaries. Namely, denoting by rj the 
elements of {0, 1}[~^'^] and by / functions on {0, 1}[~^'^], 

Lb,±fiv) ■■=^11 D^rjixMv^-^) - fiv)] (3.1) 

D+ri(x) = (1 — rj{x))rj{x + 1) ■ ■ ■ ri(N), 
D_7]{x) = ?7(x)(l - r]{x - 1)) ■ ■ ■ (1 - ?7(-A^)) 



where j > is a parameter of the model, /+ = {N — 1, A^} and /_ = {— A^, — A^ + 1} (in 
[1], [S], [6] I± consist of K sites, here we restrict to K = 2 only for notational simplicity). 
Thus Lb,+ adds a particle at rate j^ in the last empty site (if any) in /+ while at the same 
rate Lb.- takes out the first particle (if any) in J_ at same rate 1. 

Motivations for this model can be found in previous papers, [1], [5], [6], where we have 
studied the hydrodynamic behavior of the system and the profile of the stationary measure 
as A^ — 7- oo. The analysis in the above papers does not say what happens for the process 
after the hydrodynamical regime, i.e. at times longer than A^^. This is the aim of the 
current paper where we study the spectral gap, determining the time scale for reaching 
equilibrium. 

We use the same notation as in the previous section with fiNSt, t > 0, the law at time 
t of the process with generator L starting from /itv: 

fiNSt [f] = fiN [e^V] • (3.2) 

If j = i.e. L = Lq the sets {^ rj^x) = M}, < M < 2N + 1, are invariant so that the 
process is not even ergodic. However, the presence of Lf,, even if "small" due to the rate 
j/2N, changes drastically the long time behavior of the system and it is therefore crucial 
in the computation of the spectral gap. Our process, like the one in the previous section, 
is uniformly Doblin; there is therefore a unique stationary measure fi% and fl2.2p holds in 
the present context as well. We prove the analogue of Theorem [1] 

Theorem 2. There are c and b > independent of N so that 

Wf^NSt — /i/^ll < cNe' *, for all initial measures /ijv and all t > 



Theorem |2] is the main result in this paper and it will be proved in the next sections. 

The rate A^~^ in the exponent in (13. 3p cannot be improved, as can be easily seen 
by bounding from below the probability that an initially existing discrepancy does not 
disappear by the time A^^. 

The result is in several respects surprising: the spectral gap in fact scales as A^~^ just 
like in the stirring process (i.e. with j = 0) restricted to any of the invariant subspaces 
{V- "^vi^) = ^}- The result says that in a time of the same order the full process 
manages to equilibrate among all the above subsets according to /x^; also, the time for this 
to happen scales in the same way as for the process of the previous section, where however 
the birth-death events are not scaled down with A^ as in Theorem [2J 

We do not have sharp information on yU^. In |6] we have proved that the set Ai of 
all probability measures on {0, 1}[~^'^] shrinks after a time of order A^^ to a smaller set 
A^TV but we have no information on the way it further shrinks at later times. All measures 
in Ai^ are close to a product measure 'Jn, meaning that the expectation of products 
?7(xi) ■ ■ ■77(x„) are close (the accuracy increasing with A^) to those of 77V, for all n-tuples 
of distinct sites xf, n is given but it can be taken larger and larger as A^ increases. We also 
know that the expectations 77v[^(a;)] are close to p^^{x/N), where p^^{r), r G [—1, 1], is the 



stationary solution of the limit hydro dynamic equation; it is an increasing linear function 
andp"*(-l) = 1 - p'^'il) > 0. 

We thus know that fi% is close (in the above sense) to the product measure 7Ar, but 
that is all, which does not seem detailed enough to apply the usual techniques for the 
investigation of the spectral gap using equilibrium estimates. We proceed differently, our 
proof of Theorem [2] follows along the lines of the much simpler Theorem [1] It relies 
on a careful analysis of the time evolution, exploiting stochastic inequalities, as in the 
previous section. We thus consider a coupled process on X^, see ( 12. 5p . which again starts 
from ?7*^^)(x, 0) = 1 and ?7'^^'*(x,0) = for all x E [— A^, A^]. The process is defined in 
such a way that the marginal distributions of r]^^^ and ry^^^ have the law of process with 
generator L. By the definition of X^, f]^^^ > f]^"^^ at all times (order is preserved) and the 
proof of Theorem [2] follows from an estimate on the extinction time of the "discrepancy 
configuration" r/^ = r/(i)— r/(*'). We shall in fact prove that there are c and b > independent 
of A^ so that 

N 

^ P[r/^(x,t) = l]<cAre-''^~'* (3.3) 

x=-N 



4 The coupled process 

Throughout the sequel we shall use the following: 

Notation, e := A^"^ for r/ = (r/(^), r/^^)) g X^ as defined in ([23]), and x G [-A^, A^], 

?7^(x) = ri^^\x) — 7]^'^\x), rji{x) = rj^^' (x)?]^"^' (x) , ^o{x) = (1 — ?7*-"'^-'(x))(l — 7]^'^\x)), 

(4.1) 
ri^,r]i,r]o are all in {0, 1}' ^'^^ and r]^ + r]i + r]o = 1. Thus (14. ip establishes a one to one 
correspondence between X^ and {7^, 1, 0}^"^'^' and. By an abuse of notation, we shall 
denote again by 77 the elements of {7^, 1, 0}^"^'^', thinking of ?7^, ?7i, r]o as functions of rj. 
We may then say that a 7^, 1 or 0-particle is at x according to the value of 77(0;). 

Definition. Call L'q the stirring generator acting on functions on Xn (defined as in (12. ip 
with T] replaced by rj) and let Lc = L'q + 2nLi, Li = Lr + Li, be the generator acting on 
functions on X^, where Lrf is defined as 



N 



i=N-l 

N 

+ A(r/,A^)[/(r/^'i'^^°'^'^-i)-/(r/)]+ J^ B{rl,^)[f{rf''n - fiv)] (4.2) 



=N-l 



and where r/'^-'''* changes from a to 6 the value of 77 at site i if 77(2) = a, and 77°^'*'* = r/ 
otherwise, and r/^.i.^;0,^,^-i = (^^,i,JV)0,^,jv-i^ 

/^(r/,iV) = r/^(iV)[l-r/o(iV-l)], Z}(r/, iV - 1) = r/^(7V - l)r/i(iV) 

A{ri,N) = r/^(iV)r/o(iV-l) 

5(r/,iV) = r/o(iV), 5(r/,iV-l)=r/o(iV-l)r/i(iV). 

Thus Lr describes three type of events all occurring in J+: 

• D-events: a 7^-particle becomes a 1-particle. 

• A events: a 7^-particle becomes a 1-particle and simultaneously a 0-particle becomes 
a 7^-particle. 

• S-events: a 0-particle becomes a 1-particle. 

L/ is defined analogously by changing J+ into /_ and 770 with 771. One can easily check 
that 

LJ = Lg, whenever /(//) = ^(t/®), z = 1, 2, (4.3) 

L the generator in Section [31 Thus the process generated by Lc is a coupling of two 
processes both with generator L and that L preserves order. (This is just the standard 
basic coupling, as in |9]; see also Proposition 3.1 of [6].) 



5 Graphical construction 

Following the so called Harris graphical construction, we realize the coupled process in a 
probability space {Q, P) which is a product of several Poisson processes. 

Definition. The probability space {Q,P). The elements u &Q have the form 



u 



(^t(-\x e [-N,N- l];t(^'±^) t{D,±N).^iDMN-l)).t{B,±N).^(B,±iN~l))^ 



where each entry is a sequence in M+ whose elements are interpreted as times. The entries 
are independent Poisson processes: each one of the t'-^-' has intensity 1/2, all the others 
have each intensity ej'/2. 

With probability 1 all times are different from each other and there are finitely many 
events in a compact. For any such w G f2 we construct piecewise constant functions 
?7i(x, t; u), 7/o(x, t; u), r]^{x, t; u), as follows. The jump times are a subset of the events in 
the above Poisson processes, more specifically at the times t = tn we exchange the content 
of the sites x and x + 1 (i.e. we do a stirring at (x, x -|- 1)); the other jumps are: 



• At the times t = tn ' the configuration is updated only if r]^{±N,t~) = 1 and 
?7o(±(A^ — l),t~) = 1 and the new configuration has T]-t{±{N — 1),^'*') = 1 and 
?7i(±A^, t+) = 1; the values at other sites remain unchanged. 

• At the times t = tn then the configuration is updated only if r]^{±N,t~') = 1 
and rjo{±{N — l),t~) = and the new configuration has rii[±N,t'^) = 1; the values 
at other sites unchanged. 

• At the times t = tn the configuration is updated only if ?7^(±(A^ — 1), t~) = 1 
and r7i(±iV, t~) = 1, the new configuration having ?7i(±(A^ — l),t^) = 1; the values 
at other sites unchanged. 

• At the times t = tn ' the configuration is updated only if ?7o(±iV, t~) = 1, the new 
configuration having rii{±N,t~^) = 1; the values at other sites unchanged. 

• At the times t = tn ' ~ the configuration is updated only if r]i{N,t~) = 1 and 
?7o(±(iV — l),t~) = 1, the new configuration having ?7i(±(iV — l),t"*') = 1; the values 
at other sites unchanged. 

We take initially ?7^(x, 0) = 1 for all x, then the variables ri{x, t; u) defined as above on 
{Q, P) have the law of the coupled process defined in Section |H 

Definition. Labeling the discrepancies. By realizing the process in the space {Q, P) we 
can actually follow the discrepancies in time. Indeed consider the discrepancy initially at 
a site z G [—N, N]. Then the discrepancy will move following the marks of u. Namely it 
moves at the stirring times, i.e. it jumps from a; to x + 1 (or from a; + 1 to x) at the times 
t G t^^\ Moreover it jumps from A^ to iV — 1 at the times in t(^'^) (if r]Q{N — 1) = 1) and 
analogously from —N to —N + 1 at the times in t^^~^^ (if rii{—N + 1) = 1). Finally we 
say that the discrepancy dies (and goes to the state 0) at the times t^^''^^^ tiD,±iN-i)) (^jf 
the conditions for the event are satisfied, as explained in the previous paragraphs). 

We thus label the initial discrepancies by assigning with uniform probability a label in 
{1, .., 2N + 1} to each site in [—N, N] and call {zi, . . . , Z2n+i) the sites corresponding to 
the labels 1, . . . , 2N + 1. This is done independently of u and by an abuse of notation we 
still denote by P the joint law of u and the labeling. Since initially all sites are occupied 
by discrepancies we may interpret Zi as the position at time of the discrepancy with 
label i. In particular at time the probability that Zi = x is equal to 1/{2N + 1). Given 
a; e fi we follow the motion of the labeled discrepancies as described above and define 
accordingly the variables Zi{t,u;) which takes value in {[— A^, A^] U0}. Thus the set Z{t,u) 
of all Zi{t, cj) 7^ is equal to {x : ij^lI^x, t; u) = 1}, so that: 

P[^r]^{x,t)>0] = P [there is i : ^i(t,w) ^0] 

X 

< J2p[z,{t,u)j^(n] 

i 

= (2A^ + l)P[2i(t,u;)^0] (5.1) 
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the last equality by symmetry. 

Obviously P[zi(t,u) ^ 0] does not depend on the labels of the other ^;-particles so 
that we may and shall describe the system in terms of a random walk zt = zi{t,u) in a 
random environment r]t G {7^,0, 1|[~^'^1\^« when Zt ^ (i.e. it is alive); when Zt = ^ then 
rjt E {y^,0,l}^~^'^\ but since we want to study P[zi{t,u) 7^ 0] what happens after the 
death of z is not relevant. 

We have reduced the problem to the analysis of the extinction time of a random walk 
in a random environment: the problem looks now very similar to the one considered in [7], 
the only difference being that the environment has a more complex structure with three 
rather than two states per site. But the procedure is essentially the same as we briefly 
sketch in the sequel. 



6 The auxiliary random walk process 

Once the initial condition [z,!]*) has been fixed, we can consider an auxiliary time depen- 
dent Markov process (5^) as in [7], whose extinction time has the same law as that of the 
true process {zi(t)) of the previous section. The transition rates for zt are given by the 
conditional expectation of the transition rates of {zi{-)) conditioned on zi{t). Thus they 
depend on the law of the full process and hence on the initial datum {z,r]*). This time 
dependent generator Ct is given in (16. 3p below, and satisfies: 

E,[Ctf{zt)] = E,,,' [L(P{zi{t),T]t)] = -E,,,' [(P{z^{t),r]t)], 

where 0(2, 77) = f{z) and / : Aat U -^ M. 
Since 

+ ^l.=^Vo{N-l)[f{N-l)-fiN)] 
Li<f> = ^{l.=-;v(l - mi-N + l))[/(0) - f{-N)] + l,=_^+ir/o(-iV))[/(0) - f{-N + 1)]} 
+ ^l,=_^r/i(-iV + l)[/(-iV+l)-/(-iV)] 



we set 

d{N,t) = A.E,,,,4l-r]o{N-l,t)\zt = N], 

d{N-l,t) = J-E,,,,4v^{N,t)\z, = N-l] 

di-N,t) = J^E,,,,,[{l-vi{-N + l,t))\zt = -N], 

d{-N + l,t) = J^E,,,,,[r^o{-N,t)\zt = -N + l] (6.1) 

ciiN,t) = ^E,,,,4r]o{N-l,t)\zt = N], 

a{-N,t) = J-E,^,,,,M-N + l,t)\z, = -N] (6.2) 

and d{z, t) = Oif|z|<A^ — 1. Thus for t > we have 

CJ{z) = C'f{z) + diz,t)[fm-f{z)] + l,=Na{N,t)[f{N-l)-fiN)] 

+ U=.r,a{-N,t)[f{-N+l)-f{-N)] (6.3) 

The process Zt is a simple random walk with extra jumps from A^ to A^ — 1 and — A^ to 
— A^ + 1 with time- dependent intensity a{±N,t)] moreover it has death rate d{z,t) (rate 
to go to 0). Observe that 

d{z,t) > J^E,^^,'[ViiN -l,t)\zt = N]l,=^, 

and the analysis becomes very similar to the case treated in [7|. From the same argument 
leading to Theorem 1 therein, we have that for any initial configuration rj* and zq: 

P[ziit) 7^ 0] < ce-"'''"', 

which completes the proof. D 
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